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ABSTRACT 

Observational studies of galaxy isophotal shapes have shown that galaxy orientations are anisotropic: 
a galaxy's long axis tends to be oriented toward the center of its host. This radial alignment is seen 
across a wide range of scales, from galaxies in massive clusters to small Milky Way type satellite sys- 
tems. Recently, this effect has also been detected in dark matter simulations of cosmological structure, 
but the degree of alignment of dark matter substructures in these studies is significantly stronger than 
seen in observations. In this paper we attempt to reconcile these two results by performing high- 
resolution numerical experiments on N-body multi-component models of triaxial galaxies orbiting in 
an external analytical potential. The large number of particles employed allows us to probe deep into 
the inner structure of the galaxy: we show that the discrepancy between observed galaxies and simu- 
lated dark matter halos is a natural consequence of induced radial shape twisting in the galaxy by the 
external potential. The degree of twisting depends strongly on the orbital phase and eccentricity of 
the satellite, and it can, under certain conditions, be significant at radii smaller than the dark matter 
scale radius. Such internal misalignments will have important consequences, both for the dynamical 
evolution of the galaxy itself, and for mass modeling of galaxies in clustered environments. 

Subject headings: galaxies: clusters: general, galaxies: kinematics and dynamics, galaxies: evolution, 
methods: N-body simulations 



1. INTRODUCTION 

In our current picture of the Universe, structure evolves 
hierarchically from an initial density field that is very 
nearly homogenous. Tiny perturbations in this field 
seed the growth of overdensities, which then accrete 
more matter through gravitational collapse. As matter 
is pulled into these regions it is continually torqued by 
the surrounding large scale tidal field; these primordial 
torques are believed to be responsible for the orie ntations 
and a ngular momentum of galaxies at formation ([Peebles! 
fl969h . 

Recent work seems to indicate that tidal torques are 
an important effect at late times as well, acting on galax- 
ies as they fall into larger potential wells and merge hi- 
erarchically to form galaxy groups and clusters. These 
torques should affect a galaxy's orientation, and in fact 
recent observational studies have identified a significant 
correlation between a galaxy's position angle and its sur- 
rounding tidal field. These group and cluster samples, se- 
lected mainly from the Sloan Digital Sky Survey (SDSS) 
(|York et al]l200nh . show a distribution of galaxy orien- 
tations that is not isotropic, as had been long assumed; 
instead, the galaxies' long axes appear to preferentially 
point toward the center of their ho s ts dPereira fc Kulm 
20051: lAgustsson fc Brainerdl 120061 : iFaltenbacher et al 



2007). 



Subsequent numeric al studies have been able to re- 
produce these results (iKuhlen et afl 2007; Pe reira et al.1 
120081 : IFaltenbacher et al.1 l2008f ) for dark matter subha- 
los and to further establish that the anisotropy is caused 
by a dynamical process occurring inside the host clus- 
ter and which is strongly correlated with a subhalo's or- 

1 Columbia University, Department of Astronomy, New York, 
NY, 10025, USA 

2 Steward Observatory, Tucson, Arizona, 85716, USA 
Electronic address: pereira@astro.columbia.edu 



bital phase (jPereira et al.l I2008L herafter PBG08). This 
rules out a purely primordial origin, although some ra- 
dial alignment is seen even at large distances (« 3r„i r ) 
from the cluster; halos appear to be weakly aligned with 
the filaments that feed into cluster nodes. 

Under closer inspection, a significant discrepancy be- 
tween the observational and theoretical results emerges 
that has not yet been resolved. The alignment tendency 
in observed galaxies is much weaker than that seen for 
the halos in cosmological simulations (PBG08). While 
it is true that obs erved galaxy shapes can be noisy (see 
ISiverd et al.1 120091 e.g.), measurement error alone can- 
not account for such a large difference in alignment. If 
this discrepancy is confirmed, then one must conclude 
that the cluster potential is capable of inducing signifi- 
cant misalignments between a galaxy's own dark matter 
halo and its more centrally concentrated, observable, lu- 
minous component. Such internal misalignments could 
have important consequences for galaxy mass and lens 
modeling, and this possibility should be studied further. 

More generally, torquing from an external tidal field 
can affect the internal dynamics of galaxies significantly, 
and it could be disruptive enough to drive morphological 
evolution of galaxies in clusters, alongside other gravi- 
tational interactions such as harassment and stripping. 
Moreover, the alignments that this torquing induces be- 
tween observed galaxy shapes and the surrounding tidal 
field are a serious issue for weak lensing and cosmic shear 
studies, and they must be characterized and quantified 
more precisely before the next generation of cosmic shear 
surveys comes online. 

We need to study this process in more detail - how 
and when does the torquing occur? Are all components 
of the galaxy equally affected? Does the torquing induce 
radial shape twisting or do galaxies react as rigid bodies? 
Unfortunately, the low resolution of cosmological Simula- 
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tions makes the detailed dynamics of each individual halo 
impossible to study, and so we must turn to more ide- 
alized numerical experiments on single satellite galaxies. 
These can be run at higher resolution, and can include a 
dynamically distinct stellar component embedded in the 
dark matter halo. 

We begin by describing the observational param- 
eters of a typi cal galaxy in a typical cluster drawn from 
the sample of iPereira fc Kuhnl (|2005t ) (hereafter PK05), 
which we then use to build a numerical model of an el- 
liptical cluster galaxy. We describe in detail our method- 
ology (Ell), the cluster (SJI3]) and galaxy f fl2~5j) models 
used, the tidal truncation applied ( §2.6[) and the proce- 
dure used to initialize the particles in our simulations 
( £|2.7p . Halos in cosmological simulations are generally 
triaxial and we describe in section 12.81 how we include 
this intrinsic triaxiality in our model. We then briefly 
discuss the issue of shape determination ( fl2.9[l before 
presenting our results in sections [3] and |U 

As a first approach, we examine the behavior of a single 
triaxial dark matter halo orbiting an external potential 
in section [3j and compare these results with our previ- 
ous work on subhalos drawn from cosmological simula- 
tions in PBG08. We then introduce a stellar component 
in our models in section |4j We show that this observ- 
able, luminous component is also efficiently torqued by 
the cluster potential in both circular f H4.ll) and eccentric 
orbits f£ )4.2[) . although the resulting dynamics can differ 
substantially from the more extended dark matter com- 
ponent. The cluster tidal field has a significant impact 
on both the shape and internal dynamics of the satellite 
galaxy, and we investigate this in more detail in section 
[Sj We end by discussing the implications of this work 
(Q before summarizing our main conclusions in section 

m 

2. METHODS 

2.1. Methodological Approach and Code 

For the purposes of this work, we have chosen to study 
the behavior of a single galaxy on different orbits around 
a single cluster. We build these models to be represen- 
tative of a typical galaxy in a typical cluster drawn from 
the sample of PK05, as described in section 2.2. In or- 
der to study the internal dynamics of these galaxies we 
require a live, self-consistent N-body model of their stel- 
lar and dark matter distributions (described in §2.4). 
In contrast, the cluster is modeled as a static, spheri- 
cally symmetric analytic potential (§2.3), which makes 
the problem tractable numerically and also simplifies the 
geometry of the system. We do not expect dynamical 
friction between the cluster and galaxy to be important 
within these time scales and mass ratios, and the fact 
that the satellites remain on coplanar, energy conserv- 
ing orbits makes the dynamics of the system simpler to 
analyse. 

Once the initial conditions for the stellar and dark mat- 
ter particles ha ve been found ($2.7} we use GADGET 2.0 
(|Springell 12005) to evolve the models numerically. We 
have modified the code in order to adiabatically squeeze 
the live particles to form a triaxial distribution. Once 
the desired 3D shape is attained, we turn off the artifi- 
cial squeezing and let the model evolve in isolation for 
1 Gyr, at which point the shape has become stable. We 



then modify GADGET again in order to impose an ana- 
lytic potential representing the cluster tidal field, adding 
an extra acceleration component to each particle at each 
time step. The galaxy model is placed at the apocenter 
of the orbit being simulated, and the particle orbits are 
then numerically integrated for ~ 15 Gyr, or 3-5 cluster 
orbits. 

Our galaxy particles are initialized in isolation. If a 
strong external field is applied suddenly, the change in 
the potential will shock the particle distribution and lead 
to unphysical behavior. We therefore further modify 
GADGET to allow for the gradual "turn on" of the clus- 
ter potential, mimicking (crudely) the infall of galaxies 
from the field. The cluster mass is grown over a period of 
1 Gyr, and, in order to maintain the halo in the correct 
orbit through the process, we subtract at each time-step 
the acceleration induced by the presence of the cluster 
potential on the halo's center of mass from each parti- 
cle. Our results are generally insensitive to the details of 
this initialization: the overall alignment behavior of the 
halo is reproduced with or without a slow turn-on of the 
external potential. 

2.2. A Typical Galaxy in a Typical Cluster 

The observational sample described in PK05 spans a 
wide range of cluster and galaxy masses. All analyses to 
date have shown that the degree of radial alignment does 
not depend significantly on ind ividual galaxy or c luster 
properties (jPereira et all 120081: iKnebe et al.l l2008). and 
we therefore choose to build a single model representing 
a typical galaxy orbiting a typical cluster in that sam- 
ple, where we interpret typical to mean an object with 
parameters corresponding to the median of the entire 
population. 

The median redshift of the cluster sample is z = 0.1 
and the median X-ray luminosity Lx — 3.47 x 10 44 ergs/s 
(0.1 - 2.4 keV). What does the typical elliptical galaxy 
in these clusters look like? The distribution of R band 
absolute magnitudes for the galaxies in these typical clus- 
ters has a median Mr = —20.3. Their elliptical surface 
brightness profiles are generally well fit by a de Vau- 
couleu rs law over a wide range of radii (see, e.g.. lBurkertl 
(Till)): 
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where R e is the effective radius, and is roughly equiva- 
lent to the half-light radius of the galaxy. The median 
effective radius for our sample galaxies as fit by the SDSS 
data analysis pipeline is R e — 3.25 kpc. 

2.3. Cluster model 

We model the host cluster as a static external potential 
in which our galaxies orbit. We assume an NFW profile 
for the mass distribution of the cluster, and infer the 
mass and concentration of the profile from observable 
parameters of the cluster sample. 

There are many observational studies aiming to di- 
rectly relate X-ray luminosity to total mass for group 
an d cluster siz e d hal os. We use the scaling relation 
of lRvkoff etaLl pOOl who fit ROSAT X-ray luminosi- 
ties and weak lensin g masses for the MaxBCG catalog 
(|Koester et al.ll2007t ) with the following parametrization: 
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Fig. 1. — Left: The 4:1 eccentric orbit. Right: The circular orbit. Dashed blue circles represent T200 (outer) and r a (inner) for the halo 
model. 
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(2) 

This relation was calibrated for z = 0.25, and will 
evolve slightly with redshift. For a fixed mass, we as- 
sume the evolution of X-ray luminosity with redshift is 
passive and follows Lx(z) oc p c (z) 7 ^ 6 , where p c (z) is 
the c ritical density of the universe at redshift z (jKaiserl 
[I986t ). This results in a virial mass at z — 0.1 of 
-^200 = 8.3.1O 14 /i _1 M0 for the median Lx of our sam- 
ple. 

The concentrations of dark matter halos appear to be 
correlated with h alo mass in both s imulations and obser- 
vational studies. iNeto et all (|2007l ) measured this corre- 
lation for a large sam ple of z = halos dr awn from the 
Millenium simulation (|Springel et al.ll2~005f) . The sample 
spans three decades in mass and is best fit by: 



c 2 oo = 5.26 
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Halo concentrations are also observed to decrease with 
redshift, roughly along c(z) = co/(l + z), where Co is 
the concentration at z = 0. Both observed and simu- 
lated relations show significant scatter at any particular 
mass or redshift value, a consequence of the diverse ac- 
cretion histories of cosmological halos. We adjust this 
z = mass-concentration relation to the median red- 
shift of our sample and estimate a concentration for our 
typical cluster of c « 4. The parameters for the clus- 
ter NFW potential are then: M 20 o = 8.3.1O 14 fr _1 M , 
c = 4.0, r s = 532.8 kpc and r 2 oo = 2131.1 kpc. 

2.4. Characteristic Orbits 

According to cosmological simulations of large scale 
structure, most galaxies are fed into clusters through 



filaments. This leads to a radial distribution of or- 
bits in the outskirts of clusters, which progressively 
becomes more isotropic closer to the cluster center. 
Nevertheless, even for a perfectly isotropic orbital dis- 
tribution the average eccentricity of these orbits is 
high, with typical apocent er to pericenter ratios of 4:1 
(|van den Bosch et al.lll999l) . Dynamical friction progres- 
sively shrinks the orbits, but with little circularization 
- the overall measured eccentricities remain constant. 
This picture is consistent with most observational stud- 
ies, whic h find that galaxy orbits in clusters appear 
isotropic (iHwang fc Ledl2008D . although the re have been 
claims (see, e.g. iBiviano fc Poggiantil (|2009h 'l that differ- 
ent types of galaxies appear to have orbits drawn from 
different distributions, with early-types appearing gen- 
erally isotropic, and late-types biased towards a more 
radial distribution. 

The simulated clusters presented in PBG08 also show 
large orbital eccentricities in their substructure, with an 
average apocenter to pericenter ratio of 4:1, so we take 
this as the eccentricity of our representative orbit, with 
r a po = ^oo.host = 2131.1 kpc and r pori = r 20 o,host/4.0 = 
532.8 kpc . 

As a first approach, we will simplify the system even 
further and place the galaxy on a circular orbit, in order 
to first understand how our galaxy reacts to an external 
tide of constant magnitude. The radius of the circular 
orbit was chosen so that the orbit has the same energy as 
the representative eccentric orbit, i.e. r c ; r = 5r pcr i/2 = 
1332.0 kpc, and t olh » 5.2 Gyr (see figQ}. 

2.5. Elliptical Galaxy Model 

We want to create a physical model of a galaxy which 
has an absolute magnitude of M r ps —20.3 and is well fit 
by a de Vaucouleurs profile with R e w 3.25 kpc. 

The galaxy luminosity in the r band is then: L ga i,r — 
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4.76 (iBlanton et al.ll2003l) . According to the study by 
iKauffmann et al.l (12003) mass-to-light ratios vary steeply 
with absolute magnitude for galaxies observed with SDSS 
filters. For a galaxy with M r w -20.3, (M/L) r = 3.16, 
which gives a stellar mass for our galaxy of M ga i rj 
3.3 x 10 10 M Q . 

The de Vaucouleurs profile has no analytical 3D de- 
projection, but there are a few two-power density profiles 
that give projected light distributions that follow the r 1 / 4 
law at most radii. We pic k one of the mo st commonly 
used, the Hernquist profile (jHernquistll 1 9901 ) . since it has 
a simple analytical form both in density and potential 
and is finite in mass: 
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In this equation, a is the scale length, and it is related 
to the projected effective radius, R e by R e = a(l + y/2). 
Combining this with our estimate of the total stellar mass 
completely defines the stellar density profile. 

We model the dark h alo with the more t heoretically 
motivated NFW profile (jNavarro et al.lll997l) : 
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Here, p c is the cosmological critical density, S c is a char- 
acteristic dimensionless density and r s is the scale radius. 
The NFW mass profile diverges logarithmically with r, 
which means that: a) in practice the halo must be trun- 
cated at a large radius (this procedure will be described 
in £)2.6[) . and b), the mass of the halo must always be 
defined within a specific radius. 

There are a few different masses commonly quoted. 
Here, we define the mass of the halo as M200 = 200 x 
pc^itr^QQ where p c is the cosmological critical density, 
and T2oo = cr s , where c is the concentration of the dark 
matter halo. These definitions imply: 
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which means that, for a given M200, there is a single free 
parameter in the NFW profile, which can be expressed 
in terms of 5 C or c. 

Of course, the precise shape and total mass of our 
galaxy's dark component is not well constrained by our 
observations. We make a few simplifying assumptions: 
We require the circular velocity profile of the galaxy 
to be approximately flat over the observable range r ~ 
0.5 — 3.0i? e , in order to agree with the dynamical m odels 
of elliptical galaxies of, e.g. iGerhard et al.l (|2001[ ). and 
we use the mass-concentration relation of equation [3] to 
calculate c. 

These constraints limit the total mass of the halo to a 
narrow range, M200 ~ 6 x 10 11 M Q . From the mass all 
other parameters follow: Our fiducial model has M200 = 
6 x 10 n M Q , c = 7.9, r s = 22.6 kpc and r 200 = 178 kpc. 

2.6. Tidal Truncation 

The NFW profile is nothing more than a fit to halo 
density over a certain radial range. At large radii, the 



profile becomes unphysical, with the mass diverging log- 
arithmically. It therefore follows that even isolated N- 
body halos, when drawn from an NFW density profile 
must be truncated at some arbitrary large radius. The 
presence of the external tidal field in our model sets a 
natural scale for truncation, since it will strip away most 
material which lies at radii larger than the tidal radius, 
it- We approximate r f as the Roche radius of a spherical 
system under the distant tide approximation: 



/ m \ 1 / 3 
n = {ui) ^ 



(7) 



where Rq is the distance to the cluster center, m is the 
mass of the satellite and M the mass of the host enclosed 
at radius Rq. 

The exact form of the truncation is unimportant for 
our purposes, but it cannot be too sharp. Simply setting 
p = at r > r t would lead to unphysical behaviour in 
the satellite model. Instead we apply a smooth cut-off to 
the radial density profile: 
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Galaxies on different orbits will suffer different degrees 
of stripping. We choose our truncation radius so that, for 
the galaxy on the circular orbit, 80% of the galaxy's ini- 
tial mass lies within r t (r c i r ) w 80 kpc. The same galaxy 
at apocenter on an eccentric orbit will have 95% of its 
mass enclosed within r t (r apo ). These requirements set 
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2.7. Particle Initial Conditions: Eddington Inversion 

Once the density profile of stars and dark matter is 
specified, particle positions can be drawn to match these 
radial profiles using sampling techniques. However, a full 
set of initial conditions for N-body simulations requires 
the velocities of each particle to also be specified, and 
these must be such that the system is already in equilib- 
rium. 

In the case of a particle system that is described by 
a known, spherically symmetric potential, $(r), a dis- 
tribution function, /, can be constructed that depends 
solely on the two integrals of motion: the binding energy, 
E, and the angular momentum vector, L, per unit mass. 
For a non-rotating, isotropic system (which we assume 
here), L = 0, and / = f(E). The distribution func- 
tion can then be calculated using Eddington's formula 
(|Eddingtonlll916l) : 



f(E) = 



1 



V8ir 2 



E d\ 



dtp 



dp 



# 2 v^=? Ve 

where p(r) is the density profile, and ip(r) = — $(r) is 
the relative gravitational potential. The second term of 
the right-hand side vanishes for any sensible behaviour 
of tp(r) and p(r) at large r, and can be ignored. 

This formula can also be used in the case of two spher- 
ical distributions which are in equilibrium with each 
other, as is the case in our galaxy model. Then, the 
potential, VP, should be replaced by the combined poten- 
tial of the two components ^tot = ^halo + 'I'stars, 
and the integral is then solved separately for each radial 
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density profile. The distribution functions obtained will 
be for two particle distributions that are in equilibrium. 
The procedure to calculate f(E) involves first calculating 
the cumulative mass function for each component, then 
using these to calculate the total combined potential for 
the galaxy. This is gridded at regular intervals in log r 
and spline interpolation is used to obtain ^totM- This 
is then numerically integrated using an adaptive integra- 
tion routine (GSL: Galassi et al. 2006), and the values 
of f(E) are tabulated at regular intervals in \og(E). 

Once f(E) has been determined for both components, 
it is straightforward to populate the galaxy with particles 
using an acceptance-rejection algorithm (see, e.g. Press 
et al. 1996). Briefly, we start by drawing random values 
for R from a uniform distribution R G [0, -R max = 2i?2oo] 
and calculating $(-/?). Random values of V are then 
drawn from a uniform distribution V G [0, V max = 

y/— 2$(0)] and the binding energy, E, of the particle is 
calculated. If E > 0, the particle is unbound and re- 
jected. Otherwise, the normalized f(E) is calculated and 
compared to a random value, x, drawn from a uniform 
distribution x G [0, 1]. If it is smaller than x, the particle 
is kept, if larger, the particle is rejected, and the process 
starts over again, until the galaxy has been populated. 

The halo is populated with particles of mass m# — 
6 x 10 5 Af Q , for a total of Njj = 545420 particles, giving 
a total, truncated mass of Mh = 3.2 x 1O 11 M . The 
stellar particles have masses mg = 1.5 x 10 5 M Q . There 
is a total of N$ = 220000 stellar particles for a total mass 
of Ms = 3.3 x 10 10 M Q . This is enough particles so that 
two-body relaxation effects are unimportant for r > lkpc 
for a period T > 20 Gyrs, which is approximately the 
total amount of integration time required for each of our 
experiments. We check the stability of our models by 
evolving them numerically in isolation for 20 Gyrs and 
verifying that their shape and density profile does not 
evolve significantly. 

We use a separate particle softening length for the halo 
and bulge. For the halo, r so f t =0.1 kpc ~ (see, e.g. 

I Power et~a l. (2003) for a discussion of optima l softe ning 
length). Then, following iMcMillan fc Dehnenl (|2007ft . we 
pick a softening length for the stellar particles r so f t = 0.05 
which ensures that the maximum force exerted by every 
particle ( cx rni/r^ oit ) is the same. 

2.8. Triaxial Squeezing 

Triaxial halos are seen in most cosmological simula- 
tions and are believed to be a natural consequence of 
gravitational collapse and hierarchical merging. Unfortu- 
nately, finding conclusive observational evidence for tri- 
axial dark matter shapes will be incredibly difficult: the 
most promising current studies, which use tidal streams 
as particle tracers of the Milky Way potential, are still 
in their infancy. On the other hand, there are now many 
studies, both photometric and kinematic that seem to in- 
dicate that elli ptical galaxies have shapes that are intr in- 
sically triaxial (Lamba s et alJll992t iFranx et al.lll991l ). 

There is no simple recipe that can make a triaxial NFW 
N-body halo, since there are no known analytical distri- 
bution functions for systems that are both cuspy and 
triaxial. Some commonly used methods involve collaps- 
ing a self-gravit ating diffuse c o llectio n of particles nu- 
merically (e.g. iNorman et al.l (|1985l) ) or colliding two 



cuspy spherical N-body halos to cr eate a single cuspy 
triaxial one (e.g. iMoore et al.l ([2004)) but all these have 
the significant drawback of not being able to control the 
shape of the final galaxy with any precision. An alternate 
me thod, known as adiabati c sque ezing, was proposed 
bv iHollev- Bockelm ann et al.l (|2001| ) and refined later by 
IWidrowl (120081) . The main advantages of this method is 
that not only is it applicable to any equilibrium spherical 
distribution, one has reasonable control over the final tri- 
axiality of the system. We apply it here to our spherical 
model. 

2.8.1. Description of Method 

Following IWidrowl (|2008|) we numerically squash our 
spherical particle distribution to the desired triaxial 
shape. This is done adiabatically and respecting energy 
conservation which ensures that the system's differential 
energy distribution is preserved. The phase-space dis- 
tribution function is numerically altered, but the radi- 
ally averaged density and velocity profiles remain stable 
throughout the process. We modify GADGET to incorpo- 
rate the numerical squeezing of the live galaxy. 

In detail, at each timestep, an added acceleration is 
given to each particle which is perpendicular to its ve- 
locity, and obeys the following relations: 



_(fc(t)-fc(t))v 2 y +fc(t)v 2 z 



{fc{t)-fc{t))vl + fc{t)vl 

o 

V 



(10) 



where v is the speed of the particle and a • v = 0, as 
required. 

fc and fc determine the final shape of the halo. 
Roughly, a positive fc increases the x axis length, 
whereas a positive fc increases the y axis length. Most 
relevant for a triaxial system, fc — fc sets the difference 
between the minor to intermediate axis. In order not to 
shock the system with a sudden "turn on" at t = 0, the 
/3s increase from zero to their maximum values over a 
period Tq with a time-dependence given by 

fc = fc,max(t)(i(t/T G ) 2 -2(t/T G ) 3 ) . (11) 

fc and fc then remain constant for a time Tc before 
decreasing to zero over a time Td following the same 
relation. 

Tq + Tc + To must be long enough so that the process 
is adiabatic. In our case we pick Tq = Td = 0.5 Gyr, 
T c = 2.5 Gyr, fc = 0.48 and fc = 0.02, and we evolve 
the galaxy in isolation for another 15 Gyrs, to make sure 
that the shape is stable over the timescale relevant to our 
experiments. 

The minor axis length does not evolve significantly 
once the squeezing process is completed, staying fairly 
constant over the timescale of the simulations, apart from 
a small tendency for circularization that is due to par- 
ticle noise in the simulation. The intermediate axis, on 
the other hand, particularly for the dark matter compo- 
nent, shows significant evolution as soon as the squeezing 
is turned off, coming to a stable state a couple of Gyrs 
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Fig. 2. — Isophotal fitting procedure for a typical galaxy snapshot. From Left to Right: a) Projected surface brightness profile along a 
specific azimuth. Linear interpolation is used to calculate T25- b) Elliptical fit to r^h vs azimuthal angle, c) Elliptical parametrization is 
generally a good fit to the isophotal contour (in orange). 



later. It then evolves slightly towards larger values, in 
the same manner as the minor axis. 

2.8.2. Picking Axis Ratios 

Now that we have the ability to create self-consistent 
triaxial shapes, we must pick appropriate axis ratios for 
our model. We choose to squash stellar and luminous 
particles simultaneously, so the ellipticity of the halo will 
be related to the ellipticity of the stellar component. 

Elliptical galaxy shapes have a w i de dis tribution of pro- 
jected axis ratios. lLambas et al.l (|1992f ) showed that a 
combination of axisymmetric (purely prolate or purely 
oblate) systems alone cannot explain this distribution, 
but that intrinsic triaxial shapes are required. Their 
best estimate for the distribution from which these tri- 
axial shapes are drawn is a two dimensional Gaussian 
with b = 0.95, c = 0.55, a b = 0.35 and a c = 0.2. This 
corresponds to mean galaxy shapes of If)/ a) ps 0.8 and 
(c/a) sa 0.2 — 0.7 within a few times R e . 

On the other hand, we can place constraints on the 
triaxiality of dark matter halos from analyses performed 
of large cosmological simulations. In particular, a care- 
ful analysis o f the shapes o f subs tructure galaxies was 
performed by iKuhlen et al.l (|2007| ). They observe sub- 
structure halos to be less triaxial than their isolated 
counterparts, with mean axis ratios of (b/a) ps 0.83 and 
(c/a) ps 0.68 within r t , max ps 2.163 x r s ps 50 kpc for our 
model. 

Our fiducial model has stellar axis ratios of (b/a) ps 
0.82 and (c/a) pa 0.57, halo axis ratios of (b/a) ps 0.80 
and (c/a) ps 0.55 and is therefore a good representation 
of a typical galaxy according to these constraints. 

2.9. Shape Measurements 

We measure the 3D shape of our galaxy by calculating 
the reduced inertia tensor of the particle distribution: 

h X'"'^- ( 12 ) 

i i 

which weights particles equally regardless of their dis- 
tance to the center of the halo. This is then diagonalized 
to find the principal axes of the halo. The eigenvec- 
tors and eigenvalues of this reduced form of the iner- 
tia tensor give us the principal axes of the halo and a 



measure of their rela t ive le ngths (b/a, c/a). We follow 
iDubinski fc Carlberd (|1991[ ) and introduce two modifi- 
cations to the algorithm. First, the axes are found itera- 
tively: at each step the halo is rotated into the principal 
axis frame, the particle distribution and the inertia ten- 
sor recalculated. Secondly, instead of a spherical cut-off, 
an elliptical cropping is used proportional to the eigen- 
values calculated in the previous iteration. This returns 
eigenvalues that are a bit smaller, since it avoids biasing 
the measurement towards spherical values by imposing 
an artificial spherical cut-off, and the procedure gener- 
ally converges rather quickly. 

Once stellar particles are included in our model, we 
can start making more meaningful comparisons to the 
observed data of PK05. The position angles measured in 
that study were output by the SDSS pipeline and are the 
product of model independent isophotal fitting: For each 
galaxy image, an ellipse is fit to the isophotal contour cor- 
responding to I25 = 25 mag/arcsec 2 in the SDSS r band. 
The position angle of the galaxy then corresponds to the 
major axis of this ellipse. Most undisturbed galaxies have 
isophotes that are very nearly elliptical (with small devia- 
tions for boxy and disky galaxy types) and their contours 
are smoothly continuous so that their position angles and 
shapes do not vary significantly with radius. For these 
galaxies the radius of the isophote is relatively unimpor- 
tant. For galaxies with significant isophotal twisting, this 
becomes much more crucial, since this method only cap- 
tures the shape of the galaxy at the radius which it is 
measured. This is an important difference with respect 
to the method we have been using until now, in which we 
calculate the reduced inertia tensor of the entire particle 
distribution within a certain radius and thus measure a 
mean shape for the entire halo, unbiased by any partic- 
ular distance. 

In order to measure isophotal shapes for our simulated 
galaxy, we first project the stellar particles along the line 
of sight and use the mass-to-light ratio discussed in §2.51 
to calculate the surface brightness profile of each galaxy 
along 36 different azimuthal directions at 5° intervals. 
We use linear interpolation to find the radial distance 
at I25 for each azimuthal slice, and then fit an ellipti- 
cal profile to that contour (see figure [2]) . We find that 
our squeezed stellar distribution is generally well fit by 
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Fig. 3. — Sketch of angle definitions used in the text. The star 
represents the position of the cluster center. [Orbit not to scale!] 



ellipsoidal contours, with low residuals in the fit. 

It is important to note that our numerical experiments 
employ a static spherical potential, and that therefore all 
our orbits are coplanar, and the geometry of the problem 
is considerably simplified. In a real cluster, or a simu- 
lated cluster generated from the collapse of cosmologi- 
cal initial conditions, hosts are triaxial and dynamically 
evolving. Orbits precess and decay through dynamical 
friction, while also being observed from random angles, 
so the relevant radial angle is measured in 3d. With the 
experiments in this paper we are attempting to under- 
stand the dynamics of these substructures by watching 
them evolve in controlled and simplified set-ups. All sim- 
ulations were initialized so that the major and interme- 
diate axis of the galaxy are within the orbital plane, and 
the minor axis is normal to it. This is the most stable 
configuration, and it is maintained throughout the orbit. 
If the galaxy is initialized such that its minor axis is in 
the plane of the orbit, while its intermediate axis points 
out of the plane, these axes will quickly tend to flip, 
thereby achieving a more stable configuration (Christine 
Simpson & Kathryn Johnston, in prep). 

In order to quantify the degree of radial alignment, 
we calculate </>, the angle between the major axis of the 
galaxy (as determined by one of the 2 methods outlined 
above) and the direction to the cluster center. We also 
define the position angle (P. A.) of the galaxy as the angle 
between its major axis and the fixed x axis in the cluster 
frame, and 8 r as the angle between the line connecting 
the centers of the cluster and satellite and the same fixed 
axis x, so that <fi =P. A.—9 r (figure [3]). These angles are 
all always confined to the orbital plane, by design, so that 
the problem becomes essentially two dimensional. Our 
analysis will not directly address the issue of projection in 
real cluster observations, where each orbital plane is pro- 
jected onto the plane of the sky, because we are mostly 
concerned with the details of the physical mechanism. 
We note that, while this will dilute any attempt to mea- 
sure radial alignment in real galaxies, it cannot, on its 
own, account for the discrepancy between the observa- 
tions of PK05 and the simulations of PBG08, as shown 
in figure 8 of the latter. 

3. THE DARK SIDE: A SINGLE COMPONENT DARK 
MATTER HALO 

We begin by running a few simulations with a single 
component dark matter halo only, in order to directly 
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Fig. 4. — Radial Alignment, <j>, vs time for a triaxial halo on a 
circular orbit. <f> =0° means the major axis of the halo is point- 
ing in the direction of the cluster center, whereas 90° means it is 
tangential to it. 



compare the results from this idealized setup with the 
cosmological simulations analyzed in PBG08. The meth- 
ods used to create this model are essentially the same as 
for the dark matter component of the full galaxy model, 
except that the stellar potential is not included in the 
Eddington inversion procedure outlined in i)2.7l The 
shape and orientation of the dark matter halo in this 
section is always calculated from the reduced inertia ten- 
sor of particles (as described in $2.9)i that are within 
r ^ ''tid.cir ~ 80 kpc of the halo center, unless specified 
otherwise in the text. 

3.1. Circular Orbit 

As a first step we take the simplest geometric case and 
place the triaxial halo on a circular orbit, with the major 
axis initially pointing towards the cluster center, and the 
minor axis pointing out of the orbital plane. Figure 0] 
shows the evolution of the radial alignment angle, (f> with 
time. Strikingly, after one orbital period (T 5.2 Gyr), 
the halo is already very close to perfect tidal locking. 
The halo appears to oscillate about a stable position and 
the amplitude of these oscillatons dampens with time, 
so that by the third orbit, the halo is almost perfectly 
aligned. 

This is shown more effectively in figure[5j where we plot 
both the position angle of the halo and the angle towards 
the cluster center, 9 r . At first, the halo is torqued toward 
the cluster center, but by the time the halo has reached 
the first quarter of the orbit, at t w 1.3 Gyr, it has not 
rotated far enough to match the angular rotation of the 
orbit, and so it is torqued back to its initial PA. ss 0° 
alignment. On the second half of the orbit, this does not 
happen and the halo is continually torqued in the same 
direction. This sets up a uniform rotation throughtout 
the rest of the halo's evolution, where the halo shape 
is figure rotating at the same angular frequency of the 
orbital motion. 

A visual representation of the halo's alignment is given 
in figure [H] for the first orbit, where the swinging for the 
first half orbit is evident. There is also a significant 
amount of isophotal twisting apparent in the contours. 
This is most in evidence in the outermost contour, which 
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Fig. 5. — Position angle (black circles) and radial direction, 
r , (gray straight lines) vs time for a triaxial halo on a circular 
orbit. After the first orbit (t = 5.2 Gyr), the halo undergoes figure 
rotation at the same rate as its orbital angular motion so that the 
two vectors arc constantly aligned. 



represents a very loosely bound component, but it is also 
seen closer to the center, particularly for the first half 
orbit. The trends of figures 0] and [5] do not depend sen- 
sitively on where we place the outermost boundary for 
the inertia tensor measurement, but we will analyze this 
issue more carefully in ^3.31 We note that the galaxy 
begins with about 20% of the particles outside the tidal 
radius and loses mass as t~ - 2 , with about 50% of the 
mass becoming unb ound after 20 Gyr, consistent with 
previous work (e.g., lHavashi et aDl2003f ). 
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Fig. 6. — Projected particle density contours onto orbital plane 
for first passage on a circular orbit. Contours are logarithmically 
spaced at a constant projected density. Dashed lines indicate di- 
rection to cluster center, and thick lines indicate the major axis of 
the h alo as determined by the shape fitting procedure outlined in 
£12,91 The sequence starts with the snapshot at the far right, and 
proceeds counter clockwise, with each snapshot 0.65 Gyrs apart. 
[Orbit is not to scale!] 



3.2. Eccentric Orbit 

On an eccentric orbit, the magnitude of the tidal field 
will vary significantly. For an NFW density profile, the 
tidal force, F t increases monotonically with decreasing 
radius. At large radii, r >> r s , the increase is drastic 
(F t oc r~ 3 ). As the central regions of the halo are ap- 
proached, the tidal force continues to increase, but more 
slowly: at r ~ r Sl F t oc r _1 . Further in, within r < 0.1r s , 
F t asymptotes to a constant value. 

At the same time, the halo's orbital speed is no longer 
constant. Halos will spend most of their time close to 
apocenter, and pass through pericenter very quickly. The 
more eccentric the orbit, the more extreme the variation 
in orbital speed. The angular velocity is also not con- 
stant, so that the tidal torque from the host will vary 
azimuthally very rapidly through pericenter, but much 
more slowly through the rest of the orbit. How will this 
affect a triaxial halo's orientation? 

Figure [Jj shows the evolution of radial alignment for 
four orbital periods. The halo is mostly aligned through- 
out its orbit, except at pericenter (the vertical dashed 
lines in figure [Jj where there is strong misalignment with 
the radial direction. Immediately after pericenter the ra- 
dial angle decreases again, and swings back toward zero, 
overshoots, before coming back to settle at a low offset 
for the rest of the orbit. This behaviour is replicated al- 
most exactly for each orbit, and the halo shows no sign 
of asymptoting to a different alignment state, in contrast 
to the behaviour on the circular orbit, where after one 
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Fig. 7. — Radial Alignment, </>, vs time for a triaxial halo on an 
eccentric orbit. Dotted vertical lines represent pericenter passages. 



orbital period, near perfect locking was achieved. 

The misalignment at pericenter occurs very suddenly, 
and it is somewhat surprising how quickly alignment is 
reinstated. This trend was also observed for the average 
stacked population of halos in PBG08. However, our 
resolution was not high enough to observe the detailed 
dynamics for any individual halo. Naively, one would 
expect the effect of torquing to be largest at pericenter, 
since this is where the tidal field is strongest, but this 
overlooks the effect of the variable orbital speed, which 
must also be taken into account. 
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In PBG08 we suggested that the sudden misalignment 
could be due to the very high orbital speeds during each 
pericentric passage. The halos are moving too fast to be 
efficiently torqued, and therefore the radial angle, <f>, in- 
creases. Soon after pericenter the alignment is recovered, 
and the halo is torqued back onto the radial direction, in 
effect swinging around back onto itself. These high reso- 
lution experiments have confirmed this, and allow us to 
actually see this effect in action. Figure [5] shows the pro- 
jected particle density contours onto the orbital plane for 
the second pericentric passage on a 4:1 orbit. The halo 
travels in a counter-clockwise direction from right to left 
and we show here 5 snapshots separated by 0.3 Gyrs. 

3.3. Internal Structure and Shape Twisting 

Figures [6] and [8] show some amount of twisting in their 
isodensity contours. Certainly at larger radii, beyond 
the tidal radius of the satellite, one expects the cluster's 
tidal field to distort the shapes of the more loosely bound 
particles, eventually leading to the formation of tails and 
streamers. What happens closer to the center of the 
satellite? Is the internal triaxial structure of the halo 
perturbed by the tidal torquing, or does it maintain its 
shape through its orbit? 

We explore this in figure [9l which shows the position 
angle of the satellite as a function of satellite radius. 
More precisely, the orientation of the satellite is mea- 
sured at different radii by fitting ellipses (as described in 
§2.9 ) to logarithmically spaced projected mass contours, 
with the corresponding semi-major axes plotted in the 
figure. The position angle profile is shown for the satel- 
lite at the two orbital extremes, apocenter (pluses) and 
pericenter (dots). At apocenter there is almost no twist- 
ing observed (A P. A. < 15°), while, at pericenter, the 
shape twisting is extreme (A P.A. > 100°), even though 
the position angle varies smoothly with radius, with no 
obvious break within this radial range. 

Remarkably, this shape twisting is seen even at small 
radii, well within the scale radius of the halo (r s « 23 
kpc). It seems likely that the compact stellar component 
of the galaxy will react to the external torque in much 
the same way as the inner particles of the dark matter 
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Fig. 9. — Position angle vs a, the isophotal semi-major axis. 
Pluses are for the halo at apocenter (T = 3.9 Gyr), and dots are 
for the halo immediately after pericenter (T = 5.9 Gyr). While at 
apocenter there is almost no twisting in the projected mass profile 
A(f> < 15°, at pericenter the twisting is extreme, A P.A.> 100°. 



halo. The stars will therefore be misaligned with the 
outer halo, particularly at pericenter, and this must play 
a part in the discrepancy between the alignment seen in 
observations and that of simulated cosmological halos. 
We will explore this further in section 4. 

4. THE BRIGHT SIDE: ADDING A LUMINOUS 
STELLAR COMPONENT 

The large number of particles in these simulations has 
allowed us to delve deeper into the inner regions of the 
dark matter halo, and we have confirmed that, for a halo 
which is initially triaxial, the whole halo undergoes ap- 
proximately rigid body rotation through most of its orbit, 
with the exception of pericenter passages which induce 
significant shape twisting. However, the PK05 observa- 
tions are made on much smaller scales, within 2 or 3 
luminous effective radii. At this numerical resolution, 
the number of dark matter particles within this radius 
may not be enough to adequately sample the dynamics 
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Fig. 10. — Top: Position angle for triaxial galaxy on a circular 
orbit as a function of orbital time. The dark matter component 
(dark filled circles) responds more quickly to the tidal torque than 
the stellar component (open light circles) although after 1.5 orbits 
their behavior is indistinguishable. Bottom: The misalignment be- 
tween the dark matter and stellar orientation, AP.A. as a function 
of orbital time. 
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Fig. 11. — Position angle for triaxial galaxy on a circular orbit as 
a function of orbital time, with initial tangential alignment. Sym- 
bols as in the previous figure. If the galaxy is initially tangentially 
aligned to the cluster center, tidal locking takes a longer time to 
set in, but after 1.5 orbits, or approximately 7.5 Gyrs, both com- 
ponents are tidally locked in radial alignment. The dark matter 
component spins up much faster than the stellar particles: While 
the dark matter rotates a full 360° in the first orbit, the stars are 
slower to react to the external torque, and have only rotated by Ri 
180°. 



of this region. Moreover, dynamical models of ellipti- 
cal galaxies suggest that the phase space distribution of 
stars differs from that of the dark matter halo. A more 
complex model which includes a centrally concentrated 
stellar distribution is required, and this will also allow 
us to make more direct comparisons between our model 
galaxies and their observed shapes. 

Once stellar particles are included in our model, we 
can start simulating the observing process itself. In that 
spirit, instead of measuring average 3D shapes for our 
galaxy we will use isophotal contours to measure galaxy 
shapes at the same physical radius as the SDSS galaxies 
observed in PK05 ( -Z25 = 25 mag/arcsec 2 , corresponding 
to r 2 5 ~ 10 kpc, see ^2.9l for a description of the method). 
In addition, we choose to measure an equivalent isoden- 
sity contour for the dark matter, since this makes the 
comparison between the two components easier, and is 
also physically motivated by the possibility of observing 
projected mass maps through indirect methods such as 
lensing techniques. The isodensity contour for the dark 
matter is measured at /dm = 10.7 x 10 6 M Q kpc -2 , which 
results in an average contour radius < tdm >~ 40 kpc. 
The orientations calculated using this method are close 
to those obtained using the 3D inertia tensor, although 
not exactly the same, since they measure shapes at a 
specific radial distance as opposed to a mean shape for 
the entire halo. 

4.1. Circular Orbit 

We saw in the previous section that a triaxial dark 
matter halo on a circular orbit around a cluster poten- 
tial becomes tidally locked after just one orbit. What 
happens to the stellar particles that reside at the very 
centers of these dark matter wells? The stellar particles 
have a different phase-space distribution and may there- 
fore react differently to the same external tidal field. On 
the other hand the stars and dark matter are not com- 



pletely decoupled: dynamical friction acts between the 
two components and will work to bring the two compo- 
nents into alignment ((Nelson fc Tremainelll995l ). 

Using exactly the same orbital setup, we repeat the ex- 
periments of the previous section, but now with the full 
galaxy model. Figure [10] shows the measured position 
angles of both the stellar and dark matter components 
of our model galaxy as a function of time in orbit. The 
first thing to note is that after 6 Gyrs the behavior of the 
two components is indistinguishable: they both show fig- 
ure rotation in the orbital plane that matches the orbital 
period exactly. There is a significant difference between 
the two in the first orbital passage however, with the 
dark matter halo responding much more quickly to the 
external applied torque. The lag between stellar align- 
ment and dark matter alignment can be almost a full 
orbital period. Since massive clusters have recent forma- 
tion times, many of their members will have fallen in less 
than 5 Gyrs ago, and this will be an important dilution 
to the radial alignment signal. This lag in alignment will 
be explored in more detail in the coming sections. 

It is important to note that, if measured at the same 
radius, e.g. r 2 5, the orientations of the dark matter par- 
ticles and stellar particles are very similar, and always 
within < 15°, i.e. the induced misalignment between 
the dark matter and stellar alignment is mostly driven 
by the fact that the two are measured at different radii 
(each radius being representative of the spatial extent 
of the component). A misalignment between the stel- 
lar component at r 2 5 and the dark matter component at 
Tdm is also a misalignment between the inner and outer 
dark matter halo. The orientation of the dark matter 
halo appears to change smoothly between the inner and 
outer regions, as is evident in the shape twisting shown 
in figure |9l 

4.1.1. Dependence on Initial Galaxy Alignment 
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Fig. 12. — Radial angle for both stellar and dark matter com- 
ponent, symbols as before. The two components exhibit similar 
behavior throughout most of the orbit, except at pericentric pas- 
sages, where they respond very differently to the sudden changing 
torque. 



All our simulations so far have been initialized with the 
galaxy's major axis along the radial direction, the inter- 
mediate axis along the orbital direction, and the minor 
axis out of the orbital plane. This is not an unreason- 
able assumption, since, as we found in PBG08, galaxies 
are already preferentially radially aligned when they fall 
into clusters along filaments. Nevertheless, in a realistic 
scenario the initial alignment is not perfect and so it is 
important to understand how sensitive our results are to 
this initial orientation. 

We take here the extreme case and start one of our 
galaxies with the major and intermediate axes swapped, 
so that the galaxy's position angle is now tangential to 
the radial direction. The minor axis is still normal to the 
orbital plane. Both stellar and dark matter components 
initially rotate clockwise, as a response to the tidal torque 
(figure ITT)) . However, for both components, this motion 
comes to a stop and is soon reversed. The dark matter 
aligns more quickly than the stars, with the dark matter 
figure rotation accelerating rapidly to match the orbital 
period within 2 Gyrs. The stellar component takes longer 
to react to the applied torque - its position angle initially 
rotates more slowly in the first few Gyrs, although it 
eventually accelerates to match the orbital motion within 
7 Gyrs or 1.5 orbits. 

4.2. Eccentric Orbits 

It is interesting to consider what happens to lumi- 
nous galaxies on eccentric orbits, which are more typical 
of the cluster environment. As we saw with the single 
component dark matter halo, on an eccentric orbit, per- 
manent perfect radial alignment is unachievable, since 
the torque direction and strength vary too quickly along 
the orbit. The dark matter halo tends to respond quite 
quickly to the varying torque, except right at pericen- 
ter, where there is significant misalignment, leading to a 
strong swing back immediately after. 

What happens to the stars in this galaxy? Figure [12] 
shows us that, again, their behavior is very similar in all 
parts of the orbit except at pericenter, where, while the 
dark matter gets torqued back towards the cluster, the 
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Fig. 13. — Top: Stellar and dark matter position angles vs 
time for an eccentric orbit. Gray lines indicate the direction to the 
cluster center. While the dark matter particles appear to follow the 
orbital motion closely, varying their rotation speed to match the 
orbital velocity changes, the stellar particles appear to be locked in 
uniform figure rotation. Bottom: The misalignment between the 
dark matter and stellar orientation, AP.A. as a function of orbital 
time. In contrast with the circular orbit case, the misalignment 
here does not asymptote to zero, but is instead correlated with 
the orbital phase: |AP.A| peaks immediately after each pericenter 
passage, and remains small for the rest of the orbit, indicating 
internal alignment of the galaxy. 



stars respond very differently during and immediately 
after pericenter. 

This can also be clearly seen when we plot the po- 
sition angle of the isophotal contour vs time in figure 
[T3l Whereas the dark matter isophote follows the orbital 
motion rather closely, the stellar component essentially 
figure rotates at a steady rate, with the position angle 
changing linearly with time. This rate is maintained even 
through pericenter (especially in the first passage), when 
the tidal field is changing very rapidly. 

Is this a consequence of the specific geometry of this 4:1 
orbit, or a generic feature of satellites orbiting in NFW 
potentials? Figures [T5] and [T5] show the outcome of this 
analysis for three different orbits, each with the same 
apocenter, but apocenter to pericenter ratios of 6:1, 4:1 
and 2:1 (see the orbits in figure [T4"|) . It is evident from 
these plots that the 4:1 case was not unique. Certainly 
for orbits which are even more eccentric, the figure rota- 
tion induced in the stellar particles is even more uniform, 
and radial alignment more frequent. 

This figure rotation frequency, f2 ga i, appears to be the 
lowest rotational frequency that minimizes the energy of 
the system (by maximizing the mean radial alignment), 
and it can be estimated from the difference between the 
azimuthal (Q^,) and radial (f2 r ) frequencies of the cluster 
orbit: 

i\al = — • (13) 

This behavior is reminiscent of a spin-orbit resonant 
interaction, in that the frequencies of the orbits and spins 
satisfy a commensurability condition. The frequencies 
derived from this equation are used to produce tumbling 
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tracks which are shown on the right panel of figure [15] 
as dashed lines. For the top two panels, the tumbling 
slopes are a remarkable fit. The rotation of the galaxy's 
position angle is smooth and uniform, as if the galaxy is 
not aware of the abrupt change in the tidal field arising 
from the high eccentricity of the orbit. 

In this regime, the tumbling frequency, r2 ga i, is close 
to the angular velocity at apocenter, 4f apo , but it 
is much, much smaller than the angular velocity at 
pericenter,^ per .. The stars cannot react to the rapid 
variation in the tidal field at pericenter, responding in- 
stead to the mean torque throughout the orbit. 

As the eccentricity decreases, the mean azimuthal 
frequency remains approximately constant, but the 
variation in -$r along the orbit becomes smaller: 
ttt V^rr ~ 4, an order of magnitude smaller than 

at peri' at apo ' ° 

in the 6:1 orbit. fL a i is now much closer to -£ ■ and 

at peri 

we see that indeed the pericenter passages induce sig- 
nificant oscillations in the position angle of the stellar 
component (bottom panel of figure [Po]) . This leads to an 
overall offset from radial alignment, which is immediately 
apparent in the histograms of figure [TBI 

It seems, then, that orbits with intermediate eccentric- 
ities (r apo /r pci i ~ 2) can lead to significant stellar mis- 
alignment throughout the orbit. The dark matter halo 
is also disrupted by pericenter passages and less radially 
aligned on average throughout its orbit than in the more 
eccentric cases, but much less so than the stars. What is 
most striking in the 2:1 orbit is the overall misalignment 
between the stellar and dark matter particles. Simula- 
tions with apocenter to pericenter ratios of 3:1 and 3:2 



bridge the results shown in figure [T6j As the eccentricity 
of the orbit decreases, < \(j>\ > increases to higher values 
(until the orbit becomes almost circular). 

4.3. A Comparison with Observations 

These results indicate that when galaxies on orbits 
with varying eccentricities and different accretion times 
are combined in an observational sample of cluster galax- 
ies, such as that analysed in PK05, the average radial 
alignment measured for the stellar component will be 
much lower than that inferred for the dark matter from 
cosmological simulations. This is a direct consequence 
of the shape twisting observed in figure [TBI while on 
more eccentric orbits, the galaxies are internally aligned 
except at pericenter, for more circular orbits, this mis- 
alignment is present throughout most of the orbit. In 
the extreme case of a circular orbit (figure [TO]) , the two 
components eventually become aligned, but only after a 
couple of orbital passages - recently accreted galaxies will 
be misaligned. 

An evolving, non-spherical cluster potential will break 
the symmetry of these satellite orbits, which will affect 
the results in the previous section to some degree, but 
the general conclusions should remain unchanged. A 
more complete study that takes into account projection 
effects, evolving satellite orbits and the intrinsic noise in 
the observational measurements would be required for a 
quantitative comparison to observational results. 

5. THE IMPACT OF TIDAL TORQUING ON GALAXY 
SHAPES AND DYNAMICS 

5.1. Induced vs Intrinsic Triaxiality 
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Fig. 15. — Stellar and dark matter position angles vs time for 6:1 (top), 4:l(middle) and 2:1 (bottom) orbits. Left panels show the 
radial direction in gray lines and are folded at 90° and —90°. Righ t Pa nels show the running stellar position angle, with the dashed lines 
representing the tumbling frequency, £l ga i calculated with equation [13] 



A natural question to ask is whether the alignment 
toward the cluster center is a result of tidal stretching, 
where the alignment is induced by the tidal field, and 
the shape of the halo is elongated along the potential 
gradient, or more correctly described as a tidal torquing 
effect, where the halo, which is intrisically triaxial, can 
be treated approximately as a solid body which can be 
torqued about its minor axis while maintaining its orig- 
inal shape. We will see that in fact, neither simplifying 
description is adequate, and a better physical description 
lies somewhere between these two extremes. 

As a first test, we analyze the shape evolution of an 
initially spherical galaxy which is placed on a circular 
orbit around the cluster. In the bottom panel of figure 
\T7\ we see that while the projected axis ratio of the stel- 
lar component remains very close to spherical (b/a &s 1), 
the dark matter halo's shape is significantly altered by 
the external field. As a result of this, the dark matter 
halo shows instant radial alignment (since the elongation 
caused by the tides is always along the potential gradi- 



ent), and it maintains this induced alignment throughout 
its orbit. In contrast, the shape of the stellar particle 
distribution, being much more tighly bound, is mostly 
unaffected by the tidal field, and no radial alignment is 
induced. 

A simplifying, but useful, description is to analyze the 
system in terms of its characteristic frequencies. The 
outer particles in the galaxy orbit the galaxy at approxi- 
mately the same frequency as the galaxy orbits the clus- 
ter, since, by definition, the mean density of the satellite 
within the tidal (or Roche) radius is close to the mean 
density of the cluster within the galaxy's orbital radius 
{Pgai = Spduster)- For these particles, the timescale for 
the azimuthal variation in the tidal field is comparable 
to their own orbital frequencies, which means that their 
orbital structure will change significantly as a reaction to 
the external torque. 

At the other extreme, the stellar particles at the very 
center of the galaxy orbit at much higher frequencies - 
to these particles, the variation in the external tidal field 
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Fig. 16. — Left Panels: Stellar and dark matter radial angles vs time for 6:1 (top), 4:l(middle) and 2:1 (bottom) orbits. Right Panels: 
Histograms of relative frequencies of ||0|| for the dark matter and stellar components 



is essentially an adiabatic process. In a slowly varying 
potential, the action integral Jj = J Pidqi, where Pi are 
the conjugate momenta and qi are the generalized coor- 
dinates, is an adiabatic invariant, a conserved property 
of the particle orbit (Binney & Tremaine 1987). In this 
limit, the orbital structures that preserve the galaxy's 
triaxial shape are unaffected by the slowly varying ex- 
ternal torque, and the galaxy will react to the torque as 
a rigid body, figure rotating around its minor axis in re- 
sponse to the outside torque. Of course, most particles 
in our galaxy are in an intermediate regime, where some 
of the energy from the tidal torque will go into heating 
the particle orbits, and some into a coherent rotation of 
the triaxial shape. Evidence of this is seen in the evolu- 
tion of the galaxy's shape, which we discuss in the next 
section. It is important to note that a full particle orbital 
analysis would give us a better understanding of how in- 
dividual orbit families react to the torquing, and should 
be undertaken in future work. 

5.2. Shape Evolution in Substructure 



Figure [18] shows the temporal evolution of the triaxial 
galaxy shape on the 4:1 orbit around the cluster, and 
its correlation with orbital phase. The dark matter be- 
comes quite elongated at pericenter passages, as is ex- 
pected for particles that are loosely bound: at pericenter 
the tidal radius for the galaxy shrinks momentarily down 
to r t ~ 50 kpc. At the same time, pericenter passages 
make the stellar component rounder. This is probably 
mostly due to tidal shocking. The quick change in the 
potential around pericenter causes orbital mixing, heat- 
ing the stellar orbits, which will puff the galaxy up and 
make it more spherical. But what drives the increase 
in the stellar component axis ratio for the first 2 Gyrs? 
The shape of the galaxy is stable in isolation, and the 
cluster potential is turned on adiabatically, so the loss of 
ellipticity must be somehow associated with its orbital 
motion. 

This becomes clear when we look at the temporal evo- 
lution of axis ratios for a circular orbit. In this case, there 
is no tidal shocking, the magnitude of the tidal torque 
is constant, and its direction varies uniformly with time. 
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Fig. 18. — Evolution of axis ratios with time for the 4:1 orbit for 
dark matter (dark circles) and stellar particles (open red circles). 
At each pericenter passage (vertical dashed lines) , the dark matter 
becomes quite elongated, whereas the stellar component becomes 
more spherical. 
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Fig. 17. — A spherical galaxy on a circular orbit. From top 
to bottom: The position angle, radial angle, and projected axis 
ratios of the dark matter (filled dark circles) and stellar components 
(open light circles) as a function of time. Both components are 
initially perfectly spherical, but while the stellar component retains 
its shape, the more extended dark matter distribution is sheared by 
the tidal field and becomes significantly elongated. This elongation 
is at all times in the direction of the cluster center, so that the DM 
component is essentially in perfect radial alignment from t = 0. 
The stellar component orientation appears random, as is expected 
for a spherical distribution. 



Nevertheless, figure [19] shows similar evolution towards 
larger axis ratios in the beginning of the simulation, and 
another short period of circularization about 5 Gyrs into 
the orbit. 

By comparing this behavior with the bottom panel of 
figure [19] we can see that in fact the episodes of circular- 
ization occur only when the galaxy is significantly radi- 
ally misaligned (||r/>|| >35°). At these points the torque is 
strongest, and while it does induce figure rotation, bring- 
ing the entire galaxy back into alignment within a couple 
of Gyrs, it also will isotropize the orbits slowly and make 
the galaxy rounder, since the adiabatic approximation is 
not perfect. 

This trend has in fact been detected by 
iKuehn fc Rvdenl (|2005h in a study of galaxy shapes in 
clusters using SDSS. They find that the axis ratios of 
elliptical galaxies tend to become rounder in denser 
environments, even after correlations with intrinsic 
luminosity are accounted for. 

5.3. The importance of initial triaxiality 

How does the alignment depend on the initial triaxial- 
ity of the galaxy? A galaxy which is more triaxial should 
be more responsive to an external torque. We test this 
by creating a model galaxy that is more triaxial (by in- 
creasing Tc in $2.8\ from 2.5 to 4 Gyrs), and repeating 
the simulations of the previous section. In the case of a 
4:1 eccentric orbit (figure [20]), we see that, in fact, at each 
pericenter passage the position angle of the more triax- 
ial galaxy reacts to the tidal torque more strongly. This 
leads to some oscillation away from the uniform rate seen 



in the less triaxial case, although the mean rotation rate 
remains the same. These oscillations grow in amplitude 
with each orbital passage. 

Also notable is the fact that the process of sphericaliza- 
tion seems to be less effective overall for the more triaxial 
galaxy. The relative increase in axis ratios appears to be 
smaller (see bottom panel of figure [20]) . This fits in with 
the previous picture: both galaxies experience the same 
tidal field at pericenter, but, while the fiducial model ab- 
sorbs most of that energy into random orbital motions, 
increasing its sphericity, the more triaxial system reacts 
to the torque by figure rotating towards the cluster cen- 
ter, gaining coherent angular momentum. 

5.4. Isolated Evolution of a Tumbling Galaxy 

It seems clear by now that the tumbling motion in- 
duced in our triaxial galaxies by the external potential is 
not just an instantaneous tidal deformation, but rather a 
true torquing of the galaxy's intrinsic shape. In the case 
of a galaxy on a circular orbit (see, e.g. figure [10]) after 
a few Gyrs the entire galaxy becomes tidally locked, so 
that it tumbles over its long axis with a uniform rotation 
period. Once the galaxy becomes locked, is the cluster 
potential required to maintain the tumbling, or is the 
motion sustainable without the forcing of the external 
tidal field? 

In order to test this, we extract one of our models (the 
one which is more triaxial initially) which has been orbit- 
ing the cluster on a circular orbit for 10 Gyrs, and evolve 
it in isolation. We begin by determining which particles 
remain bound to the satellite, since though the satellite's 
initial profile was truncated, the torquing from the tidal 
field leads to continuous mass l oss throughout the orbit 
(see, e.g., lHavashi et al.l [20031 : iKampakoglou fc Benson! 
l2007h . 

The collection of bound particles is then extracted and 
evolved in isolation for another 15 Gyrs. Figure [H] im- 
mediately tells us that the stellar component maintains 
its tumbling motion: the position angle of the galaxy in- 
creases smoothly, at a frequency which is very close to 
the locking frequency for a period greater than a Hub- 
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Fig. 19. — Relation between shape evolution (top) and ra- 
dial alignment (bottom) for the stellar component of a galaxy on 
a circular orbit. Vertical dashed lines mark the transition be- 
tween regimes of large (\\<j>\\ >35°) and small radial misalingment 
(\\<t>\\ <35°). The figure suggests that periods of large radial mis- 
alignment (\\<f>\\ >35°) lead to circularization of the galaxy shape. 



Fig. 20. — Tidal torquing of two halos with different degrees 
of triaxiality. The top panel shows the position angle of the two 
galaxies overlayed on 9 r , the angle in the direction towards the 
cluster center (in gray). The bottom panel shows the projected 
axis ratios of the two galaxies on the orbital plane. Dashed lines 
signal pericenter passages. Red dots are our fiducial galaxy, and 
blue dots are the more triaxial model discussed in the text. 



ble time. The tumbling is encoded in the stellar particle 
orbits, and is now an intrinsic property of the satellite. 
There is a small sign of deceleration, since the slope of 
this line becomes progressively shallower, but this is to 
be expected, since any orbital mixing will naturally lead 
to a loss of coherent motion of the figure. 

For the dark matter component, the case is not so clear. 
Initially the dark matter particles also undergo coherent 
rotation, but after about 5 Gyrs, or one full rotation, 
the shape of the distribution becomes too spherical too 
determine shapes reliably, and the signature of coherent 
figure rotation disappears. This sphericalization is seen 
in both dark matter and stellar particles and must be at 
least partly driven by the fact that the cluster potential 
" turn-off' was done abruptly, and such a sudden change 
in potential leads to shocking, isotropizing the orbits. It 
is nevertheless remarkable that the stellar component can 
maintain such a coherent rotation over such a long period 
of time, with only a very small sign of deceleration. 

Figure rotations of dark m atter halos were measured 
bv lBailin fc Steinmetd (|2004l) from a large cosmological 
simulation. They measure pattern speeds for 288 ha- 
los and find that they follow a lognormal distribution 
centered at fi p = 0.148 h km s _1 kpc -1 with a width 
of 0.83, and ascribe the figure rotation of their halos 
to tidal torques from surrounding large scale structure. 
Their analysis excluded all subhalos, and also halos with 
significant substructure or signs of recent merging, and 
therefore applies to somewhat isolated halos that have 
not yet been accreted onto larger structures. For compar- 
ison, our galaxy in figure |2~T1 tumbles at a rate of ~ 27r/5 
Gyrs -1 w 1.7 h km s _1 kpc , which is an order of mag- 
nitude larger than their mean value. This is unsurpris- 



ing, since, on average, the tidal torques experienced by 
substructure are much stronger than those that surround 
isolated halos. 

5.5. Intrinsic Rotation 

If dark matter halos are already figure rotating slowly 
in isolation, will this make it harder for galaxies to be- 
come radially aligned once they fall inside clusters? A 
crude test of this was done by placing the galaxy de- 
scribed in the previous section, which has uniform figure 
rotation that is stable over many Gyrs, into the clus- 
ter again, but this time on the 4:1 eccentric orbit. The 
results of this experiment are shown in figure 1221 It is 
remarkable that the original rotation of the galaxy, rep- 
resented by the thick gray line in the top panel, is quickly 
decelerated by the external tidal field of the cluster, and 
that, after 5 Gyrs, the behaviour of the initially rotating 
halo is essentially indistinguishable from the models with 
no rotation. In the first 5 Gyrs, the rotating galaxy is 
considerably less radially aligned than the non-rotating 
models, and, interestingly, this leads to a higher degree 
of sphericalization (see bottom panel), as predicted in <JSJ 

6. IMPLICATIONS AND FUTURE WORK 

6.1. Dynamical Misalignments and Galaxy Modeling 

The experiments in this paper have shown that dy- 
namical internal misalignments between dark matter and 
stellar light are an inevitable result of tidal torquing in 
external potentials. The precise radial twisting profile 
will depend on the galaxy's orbital phase and eccentric- 
ity, but generic models can be constructed from these 
experiments that quantify the degree of misalignment vs 
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Fig. 21. — Triaxial galaxy, torqued in circular orbit for 10 
Gyrs, now evolving in isolation. Both dark matter (dark circles) 
and stellar particles (open red circles) tumble at a close to uniform 
frequency in the first 5 Gyrs. 



radius, and this is the subject of a future paper. 

While some shape twisting is expected at larger radii 
i r > ft,peri ~ 50 kpc) as a result of tidal stripping, the 
amount of twisting at small radii (r < r s ,dm ~ 22 kpc) is 
surprising, and will have important implications for mass 
modeling of galaxies in the cores of clusters. Strong lens 
models of clusters generally include bright galaxy mem- 
bers by adding pseudo-isothermal spherical mass distri- 
butions centered on the galaxies and truncated at som e 
radius, typically r trU nc ~ 50 kpc (|Natarajan et al.l f2009). 
A few studies have started to use elliptical projected mass 
distributions for the dark matter, under the assumption 
that halos are triaxial and aligned with the luminous 
component of the galaxy. The amount of shape twisting 
at r < 50 kpc seen in these simulations seems to invali- 
date this approach, and further work should be done to 
try and understand the impact of this twisting on strong 
lens models and cluster mass maps. 

6.2. Stacked Galaxy- Galaxy Lensing Measurements in 
Different Environments 

Quantifying the degree of internal misalignment in dif- 
ferent galaxy populations is also vital for stacked galaxy- 
galaxy lensing techniques, where internal alignment is 
invariably assumed in statist ical studies of dark mat - 
ter halo ellipticities (see, e.g., iMandelbaum et al"1l2006l) . 
One should be careful to exclude all galaxies in cluster 
environments from this type of study since this would 
dilute the measured ellipticity significantly. 

On the other hand, it may be possible to make use of 
these misalignments to make an interest ing measurement 
of dar k matter substructure in clusters. iNataraian et al.l 
(2009) were able to detect tidal truncation of dark matter 
halos in the cluster 010024+16, by stacking ACS images 
of cluster galaxies in separate bins according to their dis- 
tance from the cluster center. This type of analysis as- 
sumes galaxy dark matter halos align with the luminous 
galaxy within, and so images are rotated so that all opti- 
cal major axes are coincident before stacking. However, 
if dark matter halos are in fact more strongly aligned to- 
ward the cluster center than their luminous counterparts, 
rotating these images so that they are all oriented in the 
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Fig. 22. — Triaxial galaxy with intrinsic initial figure rotation, 
placed in eccentric (4:1) orbit around cluster potential. From top to 
bottom: Position angle, radial alignment and projected axis ratios 
of galaxy over 4 orbital periods. Gray line on top panel represents 
initial figure rotation of galaxy. The green points correspond to the 
rotating model, the red and blue are the non-rotating models from 
before. The dashed vertical lines represent pericenter passages 



direction of the cluster center before stacking will result 
in a stronger detection of ellipticity in the dark matter 
halo. This result would be an undeniable signature of 
dark matter that is impossible to explain within current 
theories of modified gravity. 

6.3. Weak Lensing Contamination 

Intrinsic alignments of galaxies are a serious contami- 
nation for large weak lensing and cosmic shear surveys, 
and one which has only recently started to receive fo- 
cussed attention from the lensing community. It was 
until recently assumed that intrinsic alignments could 
only occur between galaxies at similar redshifts, which 
can then be easily downweighted in the cross-correlation 
analysis. Unfortunately the alignment mechanism de- 
scribed in this paper has the side effect of creating align- 
ments between galaxies in widely separated reshift bins, 
since a galaxy that is radially aligned with the surround- 
ing cluster potential will have an orientation that is anti- 
correlated with the orientations of the background galax- 
ies being lensed by that same potential. This is a much 
harder contamination to account for, and requires a de- 
tailed understanding of galaxy alignments induced by 
tidal fields at all scales, not just in cluster potentials. 
Some attempts have been made recently to devise sim- 
ple models that predict galaxy align ments according to 
the s urrounding structure (see, e.g., Schneider & Bridle 
2009)), but further work is required to fully incorpo- 
rate the results and predictions of this paper into such a 
model. 



18 



6.4. Tidal Torquing of Disk Galaxies in Clusters 

Another possible consequence of the strong torquing of 
dark matter halos within hosts is disk warping. Because 
of their high angular momentum, disks will naturally re- 
sist tidal torquing more effectively than the surrounding 
dark matter halo, which will introduce a misalignment 
between the halo and disk. We have shown that tidal 
torquing affects all particles in the halos, even the most 
bound, so it is not unreasonable to expect that the inner 
shells should also feel these torques. 

There have been a few studies which attempt to ad- 
dress the i ssue of disk dynam i cs in tumbling dark mat - 
ter halos (jBureau et alj ' l 19991; iBekki k Freeman! I2002D . 
These studies find that a tumbling halo can induce sig- 
nificant warping in the disk, but they embed live disks 
in rigid, analytic halos with figure rotation, since it is 
difficult to generate N-body initial conditions for triaxial 
tumbling halos. This is probably not a a very realistic 
set-up, since, for a live halo, dynamical friction between 
the disk and halo will likely deform the halo potential, 
and bring the two into alignm ent. A recent paper by 
iDubinski k Chakrabartvl (2009) tries to overcome this is- 
sue by performing simulations on a self-consistent galaxy 
model with a live disk, bulge and halo component, which 
is then embedded in an analytic rotating potential which 
represents the tumbling outer halo. They find that this 
tumbling potential can induce significant disk warping 
and even trigger bar formation in a Milky Way like disk 
model. These studies were performed for isolated disk 
galaxies with dark matter halos that are tumbling at 
rates matching those measured in cosmological simula- 
tions. The work presented in this paper has shown that 
cluster potentials can also induce tumbling motion that 
is an order of magnitude larger in the halos of substruc- 
ture galaxies. It will be interesting to see how disruptive 
this effect will be for a cluster disk galaxy, and this is 
work that is currently in progress. 

7. CONCLUSIONS 

The internal dynamics of galaxies in clusters are af- 
fected strongly by the cluster potential field at essentially 
all radii and certainly well below the traditional tidal 
radius. It is a serious oversight to ignore tidal effects 
inside this radius when modeling the dynamical evolu- 
tion of these galaxies. Galaxy orientations are affected 
by tides in two ways: in the outskirts, tidal fields distort 
galaxy shapes along the potential gradient over a short 
timescale. Closer in, tidal torquing of stellar orbits dom- 
inates, inducing figure rotation of the intrinsic triaxial 
shape. 

In summary, our experiments show that: 

• On circular orbits, within an orbital period, galax- 
ies are aligned radially and figure rotating at the 
same rate as their orbital motion, with no signif- 
icant radial twisting of their shapes. The stellar 
component takes longer to become aligned, and 
this explains some of the discrepancy between the 
relatively weak observed signal, and the stronger 
alignment seen in the cosmological simulations of 
PBG08. 

• On eccentric orbits (r a po/?"peri > 4) the galaxy is 
mostly radially aligned throughout its orbit. The 



pericenter passage is too quick to change the stel- 
lar orientation significantly, affecting the galaxy 
mostly by heating its orbits, and driving the axis 
ratios to more spherical values. 

• On intermediate orbits (r a po/V per i ~ 2), the 
galaxy's orientation is still controlled by external 
torques such that the mean orientation along the 
orbit remains radial. However, because the peri- 
center passage is slower and more disruptive, the 
stellar component is not as well aligned. The dark 
matter reacts more quickly, and is generally radi- 
ally aligned throughout the orbit. 

• The induced radial shape twisting within the 
galaxy leads to a misalignment between the stel- 
lar and dark matter orientations, since these have 
different radial extents. This contributes to the 
discrepancy between the observed strength of the 
dark matter and stellar alignment signals in previ- 
ous works. 

• The figure rotation acquired by the orbiting galaxy 
becomes encoded in the stellar orbits, and is self 
sustaining for over 15 Gyrs in isolation. 

• Galaxies orbiting in cluster potentials tend to be- 
come more spherical with time: this is in part due 
to tidal shocking at pericenter, but it is also driven 
by tidal torquing during periods of radial misalign- 
ment. 

• Galaxies that fall into clusters with initial, intrin- 
sic figure rotation will generally take longer to be- 
come radially aligned than non-rotating galaxies, 
although this rotation will be damped by the ex- 
ternal torquing, and will eventually match the ro- 
tation acquired by the non-rotating models. 

A particle's reaction to an external time varying tidal 
field depends strongly on its own orbital frequency. In 
the outer regions of a galaxy, the orbital frequencies of 
the dark matter particles are close to the frequencies of 
the satellite's own cluster orbit, which are in turn di- 
rectly related to the variation in the external tide. In 
this regime, the orbital structure of particles is affected, 
leading to drastic changes in the shape and orientations 
of these outer layers on very short timescales. Stellar 
particles at the core of the galaxy orbit at much higher 
frequencies and thus see the cluster potential as changing 
more slowly, quasi-adiabatically. Their orbital structures 
are affected less, and their shapes are mostly preserved, 
although the orbital-averaged orientation still changes in 
response to an applied torque. 

This behavior leads to satellite shape profiles that can 
vary strongly as a function of radius. During pericen- 
ter passages, radial twists as large as 90° are observed 
within the dark matter scale radius, with correspond- 
ingly large misalignments between the stellar and dark 
matter orientations measured at different radii. These 
internal misalignments are also seen in the first few Gyrs 
of a circular orbit, before the satellite becomes tidally 
locked, as well as throughout most of the intermediate 
eccentricity orbits. 

We believe this internal twisting is the main reason 
behind the discrepancy between the radial alignment 
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seen in cosmological dark matter simulations and the 
weaker alignment observed in the SDSS studies. While 
the degree of twisting in dark matter halos will be dif- 
ficult to constrain observationally, these misalignments 
will strongly affect the mass modeling of cluster galax- 
ies, and should be studied further. 
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